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Magnetic fields are everywhere in nature and they play an important role in every astronomical 
environment which involves the formation of plasma and currents. It is natural therefore to suppose 
that magnetic fields could be present in the turbulent high temperature environment of the big bang. 
Such a primordial magnetic field (PMF) would be expected to manifest itself in the cosmic microwave 
background (CMB) temperature and polarization anisotropics, and also in the formation of large- 
scale structure. In this review we summarize the theoretical framework which we have developed 
to calculate the PMF power spectrum to high precision. Using this formulation, we summarize 
calculations of the effects of a PMF which take accurate quantitative account of the time evolution 
of the cut off scale. We review the constructed numerical program, which is without approximation, 
and an improvement over the approach used in a number of previous works for studying the effect 
of the PMF on the cosmological perturbations. We demonstrate how the PMF is an important 
cosmological physical process on small scales. We also summarize the current constraints on the 
PMF amplitude B\ and the power spectral index ub which have been deduced from the available 
CMB observational data by using our computational framework. 
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INTRODUCTION 

Many astrophysical and cosmological phenomena in the universe are affected by magnetic fields over a broad range 
of scales. Indeed, magnetic fields with a strength oi B ^ 1.0 //G have been detected [iHl] even on scales as large 
as that of galactic clusters. Such magnetic fields are frozen-in to the ionized baryons. The magnetic energy density 
scales as oc a~* while the baryon density scale as p oc a~^. Hence, we can relate the field strength to the cosmic 
density B^ (x p^. Considering that clusters of galaxies observed today would have isotropically-coUapsed relative 
to the background density, a cluster magnetic field of B ~ 1.0 pG would correspond to a primordial magnetic field 
(PMF) of order ~ 1.0 nG at the epoch of photon last scattering, z ^ 1100. 

The origin of such a cosmological primordial magnetic field has been an area of active research by many authors. 
One expects that such a field would have a random distribution of orientations and field strength. If the PMF has 
a nearly scale invariant spectrum, an origin from vector potentials generated during the inflation epoch is one of 
the best candidates [Jfll- Cosmological phase transitions could also produce a PMF with a bluer spectrum 
Several authors have also discussed magnetic fleld generation on smaller scales during or after the epoch of photon 
last scattering (z;^1100) ^[1-13]. Since each model for the generation of the PMF involves different length scales, the 
spectral index of the PMF power spectrum, correlates with the generation models. Therefore, constraining ne 
can lead to constraints on models for the generation of the PMF. 

A primordial magnetic field (PMF) affects electrons and protons through the Lorentz force. Subsequently, ionized 
baryons affect photons through Thomson scattering. In this way the PMF affects photons indirectly in the early 



universe. A number of studies in the literature |14| - |40j have analyzed effects of a fiducial PMF of order ~ 1 nG on the 
matter power spectrum as well as the temperature fluctuations and polarization anisotropies of the cosmic microwave 
background (CMB). Many others have been, also been attempting to constrain the strength of the PMF by various 
means, e.g. the non-gauusianity of the temperature fluctuations of CMbEJH^, Faraday rotation effects 44 - 47| . 
CMB anomalies i36, 48 . 49 ]. and effects on large scale structure (LSS) [sol. ISlj. These studies have have indicated 
that the PMF effects are mainly manifest on the smallest scales in the linear regime, i.e. before the formation of 
nonlinear structure. This is of interest because, there are possible discrepancies between observation and theory in 
the linear regime on small angular scales which could be attributable to a PMF. However, to accurately estimate the 
influence of not only a PMF but also non-linear effects, it is crucial to have the proper tools to constrain the PMF 
from cosmological observations. The purpose of this review is to summarize the development of such tools. 

We have developed methods ^25-29, 32] to analyze the effects of a PMF on the matter and CMB power spectra. 
The parameters characterizing the PMF have then been constrained from a flt |38] to the observational data. By 
simultaneously fltting the matter and CMB contributions, we have shown that a more comprehensive and accurate 
model for the PMF effects in cosmology can be obtained. This allows for improved constraints on the parameters 
characterizing the PMF. 

In section II of this review we will introduce the basic cosmological physics of the PMF and in sections HI and 
IV we will discuss its effects on the matter power spectrum and CMB. In section V we will discuss the results of a 
search for constraints on parameters of the PMF from observations of the matter power spectra and CMB. In section 
IV we will also explore the possibilities for new physics to be obtained from a better understanding of the PMF, e.g. 
contributions of the PMF to the gravitational wave mode (BB-mode). 

THE MODEL 

In this section we review how to model the effects of a PMF on the cosmic fluid. We assume a flat Friedmann- 
Robertson- Walker (FRW) background cosmology for the linear perturbations and adopt a conformal synchronous 
gauge as in Ref. [i^l- Hence, the line element is given [2^ by 

ds^ = a'^{T)[~dT^ + {S,j + h,j)dx'dx'l (1) 

where the x^ are spatial coordinates, a(r) is the scale factor, hij is the metric perturbation around the background 
spacetime, and t is the conformal time defined by: 

(2) 



»('' 



Here and in the following we use natural units c = h = 1 except where otherwise noted. 
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Primordial Magnetic Field 

In order to solve the perturbed Einstein Equations it is first necessary to specify tire energy momentum tensor for 
the magnetic field. We can adopt a prior upper limit on the PMF amplitude of order ~ 1.0—10 nG based upon 
observations of magnetic fields on the scale of galactic clusters as described below. Fortunately, for a PMF of order 
1.0 — 10 nG at the surface of photon last scattering, the total energy density in the PMF is smaller than the energy 
density in the temperature fluctuations of the CMB. Therefore, we can treat the energy density of the PMF as a first 
order perturbation and assume a stiff source for the time evolution. By a stiff source we mean that all back reactions 
from the fiuid onto the magnetic field can be discarded because these are second order perturbations. In this case, we 
can also assume that the conductivity of the primordial plasma is very large and that the electric field is negligible, 
i.e. E ^ 0. This "frozen- in" condition is a very good and useful approximation [l5| . 

On the largest scales the time evolution of the PMF can be decoupled from its spatial dependence, i.e., B(t, x) — 
Bo(x)/a^. This leads to the following simplified electromagnetic energy- momentum tensor, 

^™[EM,(x,r).|^ ' 

r'°[EM](x,T) =TO*=[em](x,t) =0 , (4) 
1 

When dealing with cosmological fluctuations, it is convenient to work in /c-space and denote all quantities by 
their Fourier transform convention F{k) — J cPxexp{ik ■ x)F{x), where fc is a wave number. Just as for the CMB 
temperature fluctuations, a PMF that is statistically homogeneous, isotropic and random, the fluctuation power 
spectrum can be parameterized as a power-law P{k) oc fc"^ |15l . \2m where ub is the spectral index. 

The electromagnetic stress-energy tensor in k space is given by [25| . 



-r''[EM](x,r) ^ai^—^ {2B\^)B'ix) - S^'^Bi^)'} . (5) 



T;ik)[EM]^^ I -0^!^l5]B\k')B,ik-k')- S'(fc')i?,(fc-fc')| 



(6) 



When comparing with observations, one desires a statistical measure of the fiuctuations on various scales. For this, 
it is best to work with a two-point correlation function for the PMF which then be defined by 

, (r) \nB+8 f>2 

B\k)B= (fcO) - ^ „, „+3 ^ .n„\3^ k--P^={k)5{k-U), k<kc , (7) 



where Bx = |B;!^| is the strength of the comoving mean magnetic field derived by smoothing over a Gaussian sphere 
of comoving radius A and k\ — 2tt/X (with A = 1 Mpc here). The tensor P'^^ {k) is defined by P'^^{k) = 5'^^ — k'^y /k^ . 
The damping scale from radiative viscosity provides a natural cutoff wave number kc in the PMF magnetic power 
spectrum. It is defined in Refs.js^-S]- 

We have evolved the PMF source power spectrum using the numerical methods described in Refs. [2^[2^. Using this 
method, we have quantitatively evaluated the time evolution of the cut off scale, and hence, more reliably calculated 
the effects of a PMF. 



Cosmological fluids, curvature fluctuatioins, and the PMF 

Cosmological perturbations can be either of a scalar form (fluctuations in energy density) , vector form (momentum 
fluctuations), or tensor form (T'-' fluctuations). Considering the scalar mode flrst, the ionized baryons (electrons 
and protons) in the early universe are influenced by PMF Lorentz forces. Before the epoch of photon last scattering 
(z > 1100 ), the photons are immersed in a fluid of ionized baryons along with the PMF and are indirectly affected 
by the PMF through Thomson scattering. At the same time, the matter and radiation respond to the background 
curvature fluctuations 77. We have found that it is very important to carefully account for the dependence of the 
curvature fluctuations on the conformal time r up to at least second order, ©(r^). This is because there is a cancellation 
between contributions from the primordial magnetic field and primordial radiation at first order. Fortunately, higher 
than first order 0{t) terms of 77 (r) can been included by considering the matter contributions to the scale factor 

MS, 



r + ^^^^^t' ^ ar + fir' (8) 
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where, puo and Pmo are total energy densities in radiation and matter, respectively. A model without this matter 
contributions may give some values which are mathematically close to the correct answer for a limited range of 
time and length scales. However, such models are not physically correct. In fact, by neglecting this coupling, the 
curvature perturbation 77 of the scalar mode of our previous numerical estimation [29j was too small to stabilize 
the numerical calculation for large scales and early times. This problem is illustrated in Figure [T] which shows a 
comparison (for B\ = 1.0 nG and ub = —2.9) of the CMB temperature fluctuations in the scalar mode with and 
without a consideration of the matter contributions to the scale factor. Since the previous curvature perturbation rj 
was too small and there was an instability in the numerical calculations for large scales and early times, the CMB 
temperature fluctuations for lower £ rise higher without the matter contributions to a than a calculation with matter 
contributions. In order to solve this problem, we have adopted the model of Ref. [s^ for estimating the effects of the 
PMF on fluctuations of the scalar mode in the early universe. We have also utilized adiabatic initial conditions for 
the matter contributions as in Ref. [H^. This leads to stable numerical calculations of the curvature perturbations of 
the scalar mode for all length scales and times. This is an improvement over previous numerical estimates for which 
the scalar curvature perturbations were too small to be stabilized in the numerical calculations for large scales and 
early times. Thus, our current method allows one to obtain consistent results for all scales and times of interest. 



Correlations between the PMF and the Primary Density Field 

Many authors have studied models for the origin of the PMF. However, there is little consensus yet as to the true 
origin of the PMF. Because of this, we cannot be certain of how the PMF correlates with fluctuations of the primordial 
density field. Nevertheless, one can define a parameter "s" to denote how much the power spectrum of the PMF 
correlates with fluctuations of the primordial density fleld [2^ [2^ . 

In the linear regime, the power spectra of density fluctuations in the baryons (Pb(fc)) and cold dark matter (/'cDM(fc)) 
in the presence of a PMF are as follows. 



Phik) — \^'5[b:FL](fc)'^fb:FL](''^)/ + \'^[b:PMF] (fc)^[b:PMF] (fc) 

+ 2(VFL](fc)'5[b.PMF](fc)), (9) 
-PcDM(fc) = ('5[cDM:FL](fc)'5j*cDM:FL](fc)) + ('5[CDM:PMF] (fc)'5fcDM:PMF] (^)) 

+ 2 ('5[cDM:FL](fc)^[bDM:PMF](fc)) ' (10) 

where the brackets denote the various two-point correlation functions as defined above, and da and a € ([b : FL], 
[CDM : FL]) designate the baryon and CDM density fluctuations without the PMF, respectively. While, dp and 
/3 e ([b : PMF], [CDM : PMF]) denote the baryon and CDM density fluctuations with the PMF included (without 
the primary fluctuations). In this power spectrum we normalize the cross correlation terms with the parameter s, 



'5[b:FL](fc)'5[b:PMF](fc)) = (VFL](fc)^[b:FL](^)) ('5[b:PMF](fc)<5[b:PMF](^))' (11) 

(^fCDMiFL] (^)'5[CDM:PMF] '^[CDM:FL] (fc)'^[CDM:FL] (^) ) ( '^[CDM:PMF] (^)'^[cdM:PMF] 

(fc) . 

(12) 

When 0<s<l,s = 0, or— l<s<Oin eqs. ([TT|) and (|12p . one has positive, vanishing, or negative correlations, 
respectively. 

Next, the Boltzmann equation can be utilized to flnd the equations of for the baryons as the follows. 



« , 2,2r / N , 1.2 P-fO n[EM:S](fc) ,^ 

Vb = Vb + Cgk Ob + T^aneaTyv-. - Vb) + k , (13) 

dpb apbo 47rp-yo 



where Vb and v^ are the baryon and photon velocity perturbations, Sb and 5-y are the baryon and photon density 
perturbations, Cs is the sound speed, rig is the free electron density, ctt is the Thomson scattering cross section, pbo 
and p-yo are baryon and photon densities, and H[EM:S](fc) is the square root of the power spectrum function for the 
Lorentz force (given in [23] ) . 

The Lorenz force term in Eg. ljlSp can be divided into two terms, the magnetic pressure and the tension. By 
comparing those terms, one can decide which of them is dominant in the Lorenz force. There is, however, no 
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information as to whether the magnetic pressure or the tension is dominant, and whether the direction of the forces 
from them are the same or different. Nevertheless, such information should be taken into account when it can be 
determined from a model. To see how these pieces arise, we summarize the derivation of the scalar Lorentz force term 
as given in Ref. psj . Analogous derivations of the vector and tensor two-point correlation functions can be found in 
Ref. [29[. One can first rewrite the two-point correlation function for the scalar part of the electromagnetic energy 
tensor in k space as, 

(r[EM:S](fc)rfEM:S](fc')> - A:,Flzfc™(Tj[EM](fc)r*'[EM](fc')> 

= (2^)3|n[EM:S](fc)l''5(fc-fc') , (14) 

where n[EM:S] (^) is the power spectrum of the Lorenz force that we wish to analyze. 

One can next decompose the electromagnetic stress-energy tensor in k space into two parts as follows, 

T;(fc)[EM] = 7;'(fc)[EM:l] - r;(fc)[EM:2] (15) 

T;(fe)[EM:i] = ^ / ^l5)B\k')B,ik- k') (16) 

T,'(fc)[EM.] = ^ / ^B^ik')B,ik-k'). (17) 

Correspondingly, one can define [25] r[EM:Si] = fcifc-Tj [em:1] and r[EM:S2] = fcife'''7j[EM:2]- Using Eq. ([15] - [TTl) , one can 
finally rewrite the two point correlation function for the scalar part of the Lorenz force as, 

(r(fc)[EM:S]T"*(fc')[EM:S]) = ((T'(^)[EM:S1] - T'(fc)[EM:S2])(r*(fe')[EM:Sl] " C^') [EM:S2] )) 
= (T(fc)[EM:Sl]T*(fc')[EM:Sl]) - {Tik)[-E,M:Sl]T*{k')[EM:S2]) 

- (T'(fc)[EM:S2]r*(fc')[EM:Sl]) + {T{k)[EM:S2]T*{k')[EM:S2])- (18) 

The first term on the R.H.S. can be identified as the magnetic pressure while the fourth term is the magnetic tension. 
The second and third terms are the cross correlation. A key part of the analysis of the effects of the PMF on the 
cosmic fluid is to understand the relative roles of these terms. 

Fig. [2] shows the ratio of the magnetic pressure to the tension as in the Lorentz force term |23j as a function of the 
spectral index ub- The pressure dominates when tib < —1.5, while the tension is slightly larger when ne > —1.5. 

The sign of the cross correlation terms depends upon the relative signs of the pressure and tension. In order to 
determine the relative signs of these two terms, however, one must specify a model for the generation of the PMF. 
There is as of yet no consensus model. Therefore, we decompose the factors into various possible combinations, i.e. 

S S[LF] X spF], (19) 

where 



-1, 


n < 


-1.5 


(I), 


-1, 


n > 


-1.5 


(11), 


1, 


n > 


-1.5 


(HI), 



S[LF] = <( -1, n > -1.5 (II), (20) 
and 

< S[DF] < 1 (i), 
S[DF] = (u), 

-1 < spF] < (iii). 

In the different regimes, S[lf] represents either: (I) the pressure dominated case; (II) the tension dominated case, 
where the magnetic field pressure and tension forces act in the same direction; or (III) the tension dominated case, 
where the magnetic field pressure and tension forces act in the opposite direction. On the other hand, S[df] represents 
either: (i) a positive correlation between the matter and PMF distributions; (ii) no correlation; or (iii) a negative 
correlation. Thus, if s < 0, the matter and PMF distributions could be correlated positively (spF] > 0) and the PMF 
pressure dominates in the Lorentz term (for n < —1.5). Another possibility is that the matter and PMF distributions 
negatively correlate (spF] < 0) and the PMF tension dominates in the Lorentz term (for n > —1.5) and the tension 
acts on the density field in the same direction as the magnetic field pressure. Yet another possibility is that the PMF 
tension dominates in the Lorentz term (n > —1.5), but the tension acts on the density field in the opposite direction 
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from the pressure force. In these cases the PMF effects act Uke a gas pressure to oppose the gravitational collapse 
and therefore causes the density perturbations to more slowly evolve. 

On the other hand, if s > 0, the matter and PMF distributions could positively correlate (s[df] > 0) and the PMF 
tension dominates in the Lorentz term (n > —1.5) while the tension acts on the density field in the opposite direction 
from the pressure force. Alternatively, the matter and PMF distributions could negatively correlate (s[df] < 0) and 
the PMF pressure dominate in the Lorentz term (n < —1.5). In these cases the Lorentz force from the PMF accelerates 
the gravitational collapse. After decoupling, S does not oscillate and the perturbation evolution is straightforward for 
all of the above cases. 



EFFECTS OF A PMF ON THE MATTER POWER SPECTRUM 



In this section we review the effects of the PMF on fluctuations in the matter density on cosmological scales (see 
^ for details). 



Effects of a PMF on the Matter Density Field 



Figure. [3] shows that the density fluctuations of matter are more strongly affected by a PMF for wavenumbers 
k/h > 0.1 Mpc~^|25j. This is because the energy density of the PMF, E-q, only depends upon the scale factor as 
a^^ like the photons, and unlike the photons, the magnetic field fluctuations can survive on scales below the photon 
diffusion length, i.e. below the Silk damping scale. Since the CDM and baryons interact with each other through 
gravity, the Lorenz force from the PMF can also indirectly affect the cold dark matter (CDM). This effect on the 
CDM is much smaller than that on the baryons before the epoch of photon last scattering because the density of 
baryons oscillates with that of photons and their gravitational effect on the CDM is very small. After the epoch of 
photon last scattering, the time evolution of the CDM density starts to be influenced by the baryon density through 
gravitational interaction 57[. Therefore, the Lorenz force effect of the PMF on the CDM increases with time. 

The strength of this effect is dependent upon the ratio of the baryon density to the CDM density, fib/Oc. The 
density fluctuations of the baryons are directly generated by the PMF, while the density fluctuations of the CDM are 
indirectly generated from the gravitational interaction with the baryons. Therefore, the density fluctuations of the 
CDM only grow due to gravitational instability, and their growth rate is nearly the same as the density fluctuations 
of the primordial CDM. On the other hand, the evolution of the baryon density fluctuations is not dominated by the 

gravitational potential until the gravity becomes comparable to the Lorentz force, i.e. 

n[EM:S](^: r) is the Lorentz force power spectrum [25| . This is because the baryons density fluctuations are generated 
by the PMF directly. As mentioned above, the effect of a PMF on the baryon density fluctuations is largest on the 
smallest scales (larger k values) . Hence, the density fluctuations for baryons in the presence of a PMF increase with 
the wavenumber fc 12511. 



Matter Power Spectra with a PMF 



It should be noted that the effect of a PMF on the matter power spectrum, P(k) is different from effect of a PMF 
on the matter density fluctuations S. The fluctuations of total density S can become smaller or larger depending upon 
whether the dominant effect is from the pressure or the tension of the PMF. On the other hand, the effects of a PMF 
upon the matter power spectrum are dependent on how well the spectrum of the PMF correlates with the primary 
density ffuctuations. 

Figure [3] illustrates the degeneracy between the PMF- matter-density correlation factor S[df] and the PMF spectral 
index ns- This flgure shows that a change of spectral index from ns — —2.0 to —1.0 can be offset by including a 
negative correlation. The correlation factor S[df]) however, depends upon the origin of the PMF. We will discuss the 
below the resolution of this degeneracy problem, along with the relation between the correlation factor S[Dp] and the 
origin of the PMF. 
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Negative correlation of the density perturbations 

In the case of a negative correlation between the PMF and the matter density fluctuations, the pressure of the PMF 
accelerates the evolution of the fluctuations in the matter density, while the tension of the PMF delays them. Thus, 
for n < —1.5, for which the pressure of the PMF dominates, the PMF accelerates the evolution of the matter density 
fluctuations. On the other hand, for n > —1.5, for which the tension dominates, the PMF delays the evolution. 
If we assume that the PMF would have been produced in low density regions, such negative correlations might be 
allowed. This situation, however, is difficult to realize for a causally produced PMF, because it seems natural that 
more amplitude for the PMF would be generated in regions with higher matter density and therefore higher currents. 



No correlation 



The PMF and the matter density fluctuations are uncorrelated when S[df] = 0. The PMF then increases the 
matter power spectrum independently of whether the pressure or tension dominate the PMF. If density fluctuations 



generate the PMF as in Ref. [13|, peaks in the amplitude of the PMF should lie at peaks in the pressure gradient of 



the cosmological fluids. In this case, the PMF is produced along the border between high and low density regions, i.e., 
(5 ~ 0. Eventually, such a PMF has no (or very weak) statistical correlation with fluctuations in the matter density. 



Positive correlation 



In the case of a positive correlation between the PMF and fluctuations in the matter density, the evolution of the 
matter density fluctuations is accelerated by the tension of the PMF. Therefore, the matter evolution is delayed by 
the pressure of the PMF. Thus, for n < —1.5, for which the pressure of the PMF dominates, the evolution of the 
matter density fluctuations is delayed by the PMF. On the other hand, for n > —1.5, a PMF accelerates the evolution. 
If we assume that the PMF would have been produced in higher density regions, such positive correlations might be 
allowed. 



EFFECTS OF A PMF ON THE CMB 



In this section, we review effects of a PMF on the temperature fluctuations and polarization anisotropics of the 
CMB (see [29|,|38l for details). Re garding the CMB polarization, the reader should be aware that in cosmology the 
intensity of polarized radiation is expressed in terms of two scalar fields E and B that are independent of how the 
coordinate system is oriented. E and B are the curl-free and curl-like components of the linear polarization field. 
Based upon these there are three types of fluctuations: scalar, vector and tensor. There are then four observables: the 
temperature, i?-mode, i?-mode, and the temperature cross polarization power spectra. From these, one can generate 
three power spectra TT, EE, and BB, where T is the total intensity and also three cross-spectra TE, TB, and EB. 
However, the only nonzero spectra are TT, EE, BB, and TE due to parity considerations. 

Figures |4] - [6] show that the usual (TT mode) CMB power spectra for large multipoles £ (small angular scale) are 
influenced most strongly by the PMF. The first reason for this is that the energy density of the PMF scales with the 
cosmic expansion as just like the photons. Unlike photons, however, the magnetic fields are not much affected by 
photon diffusion. As a result, the temperature fluctuations and polarization anisotropics of the CMB continue to be 
affected by a PMF even for angular scales smaller than the Silk damping scale. Another reason is that a dominant 
contribution to the CMB on smaller angular scales comes from the vector mode of the PMF. This is because, after 
horizon crossing, the scalar mode cannot evolve due to the acoustic oscillations. That is, the photons fall in and out of 
gravitational potentials but cannot grow. On the other hand, the vector mode of the PMF fluctuations can continue 
growing once inside the horizon. For higher £, the vector mode also dominates the CMB polarization from the PMF, 
and the tensor mode from the PMF diminishes [cf.figures S]-[6][i3,[2i Issl . [sit ■ The gravitational waves generated by 
the PMF can be negligibly small after horizon crossing. This is because this homogeneous solution starts to oscillate 
inside the horizon and decay rapidly [59l - l62l | . Consequently, the anisotropy spectrum is affected by gravity waves only 
on scales larger than the horizon at the epoch of photon last scattering. 

Panel (b) on figure H] and figure |6] show that the BB mode from the PMF can dominate for B\ > 2.0 nG for 
multipoles with £ > 200. On these angular scales there is also the BB mode which is converted from the EE power 
spectrum by gravitational lensing [g^ . Therefore, there is a degeneracy between the PMF and gravitational lensing 
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on small angular scales. However, since the spectrum from the effects of the gravitational lensing signal have been 
independently determined, this power can be subtracted directly. We have shown [ii 113 that the effects of a PMF 
on the temperature fluctuations of the CMB are stronger with correlations than than without, as is also the casefor 
the matter power spectrum with a PMF. This is illustrated in Figure [S] 

There is also a degeneracy between the Sunyaev-Zel'dovich effect and the PMF. Nevertheless, the effects of a PMF 
on the CMB are independent of frequency because the PMF influences the primary CMB as a background. Eventually, 
by using observations at different frequencies, it should be possible to distinguish the effects of a PMF from foreground 
effects which depend upon frequency. 



CONSTRAINTS ON THE PMF 



In this section, we review constraints on the parameters of the PMF from cosmological observations. First, we show 
that a strong constraint on the PMF amplitude and spectral index can be derived from the prior observed constraints 
on the cTg parameter. This fiducial quantity represents the root-mean-square of the matter density fluctuation in a 
comoving sphere of radius 8h~^ Mpc and as such is a measure of the growth of large scale structure. It is given by 
a weighted integral over the matter power spectrum ^63]. Second, we note that the parameters of the PMF can be 
strongly constrained by using the observational data on the CMB with the previous priors used in the analysis of 



Refs. [6J-|66[ without a PMF. In the last subsection, we will discuss the possible origins of a PMF based upon the 



deduced constraints on the parameters characterizing the PMF. 



Prior on the PMF Parameters from erg 



Because of the degeneracy of affects on the observed CMB power spectrum from various parameters, it has become 
standard in the field of CMB physics to fit the CMB parameters by a Baysian method that makes optimum use 
of the varity of independently determined constraints on cosmological parameters. For this purpose the Monte- 
Carlo Markov Chain (MCMC)j67i] analysis has become the method of choice to find the multidimensional surface of 
maximum likelihood among the many parameters. In this context we note that the PMF amplitude B\ and power law 
index ub have a strong degeneracy. Therefore, one must identify a strong prior constraint to effectively decouple these 
two parameters of the PMF. In this quest we are aided by observed constraints on the as parameter for which a value 
of 0.7 < CTs < 0.9 has been deduced from diverse observational data on linear cosmological scales [68l - [7lj |. Using this 
range as a prior in the Monte-Carlo Markov Chain (MCMC)[63 analysis, constraints on the parameters of the PMF 
can be obtained by a fit to the observed power spectra of the CMB. For some of the cosmological parameters, ilb, 
^^CDM , "-s I a-nd affect as , we also have to consider the possible degeneracy between these cosmological parameters 
and the PMF parameters. Fortunately, however, recent observations of the CMB have constrained such cosmological 
parameters on larger scales {i < 1000]l64,-ii6|. As mentioned above, the PMF mainly affects the CMB anisotropics 
on smaller scales ( £ > 1000) [H, 26, 2^. Therefore, only a small degeneracy between the PMF parameters and the 
other cosmological parameters is expected, so that fixing the other cosmological parameters at the WMAP best fit 
values is justified. 

Figure [7] and the panels in Figure [5] show the optimum PMF parameters tib and Bx for various constant values of cg 
as labeled. Since the PMF power spectrum depends upon fc^"^+'^[29[, for tt-b < —1.5, the PMF effects on the density 
fiuctuations for smaller scales decrease with lower values for riB . Therefore since trg depends upon the amplitude of the 
matter power spectrum for smaller scale, when the spectral index is near hb — —3.0, the matter power spectrum in the 
presence of a PMF is smaller for smaller scales. Thus, larger amplitudes of B\ are allowed. For nB > —1.5, however, 
the energy density power spectrum of the PMF is proportional to the cut-off scale k^^^'^ ^9^. kc is also, proportional 
to _B^i/("B+5) [^^[g|_g||_ Substituting these relations into Eqs.(15)-(17) in Ref. @, we easily obtain the following 
relation between the PMF power spectrum and the magnetic field strength, P(fc)pMF oc 

^[i4/(nB+5)]_ Considering 

Bx <^ IG, the matter power spectrum with the PMF for ub > —1.5 becomes larger for larger tib Therefore larger 
amplitudes of B\ are not allowed for larger ub- 

Since the effects of a PMF for negative and positive correlations change at ub = —1.5 (see Section p), we divide 
the discussion below of each correlation into two parts based upon whether the spectral index is greater or less than 
riB = —1.5. 
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No Correlation 

As mentioned above in Sub-Section„ in the case of no correlation between the PMF and the density fluctuations 
from primary perturbations, the effect of a PMF is to increase the matter power spectrum independently of whether 
the pressure or tension dominates the PMF. As noted above, from recent constraints on the (Ts parameter from 
cosmological observations of LSS we can exclude PMF parameters which imply dg > 1. 

Hence, we can exclude a PMF with field strength Bx ^1 nG if ub > —0.9 (Fig. [7]). We can also exclude a PMF 
with amplitude B\ ^0.1 nG if ne > 0.2. If we assume that the origin of the observed magnetic field in the clusters 
of galaxies is the PMF, we can limit the PMF amplitude to B\ ^1 nG. In this case, we can exclude a PMF spectral 
index in the range > —0.9. 



Negative Correlation 

As noted above in subsection „ in the case of a negative correlation, the tension of the PMF delays the evolution of 
the matter density fluctuations, while the pressure tends to accelerate them. Furthermore, the pressure of the PMF 
dominates for ne < —1.5, while the tension of the PMF dominates for n-e. > —1.5. Therefore, for n < —1.5 the matter 
density fluctuations are accelerated by the PMF, while, for n > —1.5, this evolution is delayed by the PMF. This 
behavior can be traced to the third terms in both Eqs. fTTI) and (fT2|). Using the allowed range of PMF parameters 
as mentioned above, we can exclude the range of PMF amplitude to B\ ^1 nG if ub > —0.81, and the range of 
PMF amplitudes to B\ ^0.1 nG if riB > 0.26(panels (a-1) and (a-2) in Figure [8]). However, as mentioned above, it 
is difflcult to model the production of a PMF in matter fields of low energy density in order to realize such negative 
correlations. 



Positive Correlation 

As noted above in Sub-Scction„ in the case of a negative correlation, the pressure of the PMF delays the evolution 
of the matter density fluctuations while the tension accelerates them. Considering the same conditions as in the 
previous section. V A 1, the PMF decreases the fluctuations of matter density for ub < —1.5. On the other hand, for 
ub > —1.5, it increases the fluctuations. In this case, we can exclude the range of PMF amplitudes with B\ ^1 nG if 
ne > —0.94, and the range of PMF amplitudes with B\ ^0.1 nG for ub > 0.13 (panels (b-1) and (b-2) in Figure (S]). 

Constraint on the PMF from the CMB 

In this section, we summarize the current constraints on parameters of the PMF from flts to the CMB and LSS 
observational data. We consider a flat (fc — 0) ACDM universe characterized by 8 parameters, i.e. {flfo/i^, ilch^,Tc, Ug, 
log{10-^^ As), At /As, \B\\, Ub}, where fih/i^ and ^Ich^ are the baryon and CDM densities in units of the critical density, 
h denotes the Hubble parameter in units of 100 km s~^Mpc~^, tc is the optical depth for Compton scattering, Us is 
the spectral index of the primordial scalar fluctuations, As is the scalar amplitude of the primordial scalar fluctuations 
and At is the scalar amplitude of the primordial tensor fluctuations. We define the tensor index for the primordial 
tensor fluctuations as nt = —(As/At)/8. For all cosmological parameters we use the same prior constraints as those 
adopted in the WMAP analysis [72] ■ As noted in the previous section, the case of no correlation is more realistic 
than the case of negative or positive correlations. Therefore, we focus here on the constraints of the parameters 
characterizing the PMF in the case of no correlation. 

Using the MCMC algorithm and available cosmological observations, we have constrained the standard cosmological 
parameters and the PMF parameters listed in Tabid. We note, however, that the MCMC algorithm tends to sparsely 
sample near a boundary, and hence find low probability there. This is why it appears that ns is constrained even in 
the limit that B\ — > 0. Of course, for B\ — there is no constraint on the spectral index at all. This shortcoming, 
however, does not negate the fact that there is a genuine minimum in the goodness of fit for a finite magnetic field and 
spectral index. In the optimum fit with the inclusion of parameters of the PMF, the minimum total changes from 
2803.4 to 2800.2 corresponding to a change in the P^r degree of freedom from 1.033 to 1.031. Therefore, a finite 
PMF slightly improves the goodness of fit even after taking into account the new degrees of freedom. The existence 
of a PMF, however, is still only of marginal significance. 
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Figures [9l and [TOl show the 68% and 95% C.L. probabihty contours in the planes of various cosmological parameters 
versus the amplitude |-Ba|, or power law index ne, along with the probability distributions. Figure. [10] and the 
bottom panels in Fig. [HI show the probability distributions for \B),\ and tib- Of particular note for this review is 
the presence of maxima in the likelihood functions for \Bx\ = 0.85 ± 1.25 nG and = — 2.371q73. Although these 
values are consistent with zero magnetic field, and thus only imply upper limits, they suggest the possibility that with 
forthcoming data (particularly for large CMB multipoles) a finite magnetic field may soon be detectable. 

These figures exhibit no degeneracy between the PMF parameters and the standard cosmological parameters. Table 
1 confirms that the standard cosmological parameters are not significantly different from those deduced directly from 
the WMAP 5yr data analysis without a PMF. The reason for this is simple. The standard cosmological parameters 
are mainly constrained by the observed CMB power spectrum for low multipoles £ < 1000 (up to the 2nd acoustic 
peak). On the other hand, the PMF dominates for £ > 1000. Hence, the PMF effect on the power spectrum is nearly 
independent of the standard cosmological parameters. 

The tensor to scalar ratio At/Ag deduced in our analysis is smaller than the upper limit At/Ag < 0.43 (95% 
CL) deduced from the WMAP 5yr data analysis without a PMF. The reason for this is that we define At as the 
tensor amplitude of the primary CMB spectrum (without a PMF). This tensor term only arises from the primordial 
gravitational-wave background produced during inflation. We combine the tensor amplitude from the PMF and At 
when we compare our tensor amplitude with the result deduced by others. The value of At by itself is comparable 
to the tensor contribution from the PMF. The value of At/As from the WMAP 5yr data is less than 0.43 (95% CL). 
Hence, our result is consistent with the previous constraints when the additional PMF contribution is included. 

The degeneracy of the PMF parameters [2^ is broken by the different effects of the PMF on both the matter 
power spectrum and the CMB power spectrum. The vector mode can dominate for higher £ in the CMB temperature 
anisotropics (29j , while the matter power spectrum becomes sensitive "2^ to the power law spectral index rte when a 
PMF is present. Additionally, the CMB fluctuations from the PMF are smaller than the primary CMB fluctuations 
for the scalar and tensor modes on large angular scales 3^ 3l|. Therefore, the tensor to scalar ratio is not affected 
by the PMF. 

Figure [TT] shows our deduced probability distributions and the la and 2(t (68% and 95% C.L.) probability contours 
for the resultant cosmological parameters, erg, Hq, Zroion, and Age. Here, Hq = lOO/i km Mpc^^ is the Hubble 
parameter, Zj-cion is the red shift at which re-ionization occurs, and Age is the presently observed age of the universe 
in Gyr. It is important to keep in mind that these parameters are not input parameters, but are output results. The 
sum of the results on Figures [9] - [Tl] provides marginal evidence that both upper and lower limits to the parameters 
of the PMF can be deduced. 

Table 1 summarizes the upper limits to the PMF parameters with input and output cosmological parameters, erg, Hq, 
Zroion, and Age. In particular we find that \Bx\< 2.10(68%CL) nC and < 2.98(95%CL}nG and nB< -1.19(68%CL) 
and < — 0.25(95%CL) at a present scale of 1 Mpc. Although previous work [2^, |2^, [2^ [s^l could obtain a less stringent 
upper limit to \Bx\ they could not constrain tib at all. Moreover, our deduced probability distributions suggests that 
a finite PMF provides the best fit. 

On angular scales smaller than that probed by the CMB the observed number density of galaxies is a better measure 
of the power spectrum. Therefore, since the PMF mainly influences the small angular scales, using the combined 
LSS observational data (2dFDR;68]), and the CMB (WMAP 5yT^, ACBARf?!, CBljli, Boomerang[7|) we can 
constrain the PMF better than in previous works which relied on the CMB data only. In particular, an upper limit 
on riB can be constrained for the first time while the lower limit to ub is approaching the la confidence level. 

The right-bottom panel of Figure [H] shows that the maximum likelihood is for a spectral index of rte = — 2.41q'7. 
It is important 2^, 77 1 to constrain rte as this parameter provides insight into models for the formation of the PMF. 



If the PMF were formed during inflation one would expect a scale-invariant value of riB = —3. The value deduced in 
Figures [51fTT] and Table 1, are thus consistent with an inflation generated PMF at the la confidence level. 

Next we discuss other generation models for the PMF during various phase transitions from the gravity waves 

produced along with the PMF. The formation of light elements during big-bang nucleosynthesis (BBN) depends upon 
a balance between the nuclear reaction rates and the expansion rate of the universe. Since gravity waves contribute to 
the total energy density they affect the expansion rate. Hence, they are constrained by a comparison between the BBN 
predictions and the observed light element abundances [1^ ■ There are also other theoretical constraints on the 
origin of the PMF. For example, in Ref. [tI] it was noted that a non-scale-invariant spectrum of the PMF is difficult 
to generate or survive if it is formed via a causal mechanism during later phase transitions. Indeed, as pointed out 



in Ref. |77l . |78| the spectral index of a magnetic fleld generated during any period of standard Friedmann expansion 



(i.e. any time except inflation) is constrained to be a positive even integer larger than two: hb ~ 2, 4, 6 This 

constraint, which applies for example to a PMF generated at a phase transition, or by second order MHD processes 
and so on, comes simply from the fact that the magnetic fleld generation is a causal process (i.e. the universe has a 
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Cosmological Parameters 


Parameter 


mean 


best fit 


nth'-' 


0.02320 ± 0.00059 


0.02295 




0.1094 ±0.0046 


0.1093 


rc 


0.087 ±0.017 


0.082 


ris 


0.977 ±0.016 


0.970 


In(lO^OA) 


3.07 ±0.036 


3.06 


At/A, 


< 0.170(68%CL),< 0.271(95%CL) 


0.0088 


|BA|(nG) 


< 2.10(68%CL),< 2.98(95%CL) 


0.85 


riB 


< -1.19(68%CL),< -0.25(95%CL) 


2.37 




PI oi9+U.Uy8 
U.01^_Q 033 


0.794 




73.3 ± 2.2 


72.8 


^reion 


10.9 ± 1.4 


10.5 


Age(Gyr) 


13.57 ±0.12 


13.62 



TABLE I. PMF parameters and ACDM model parameters and 68% confidence intervals {At/As is a 95% CL) from a fit to the 
WMAPjli ± ACBArIzI ± CBllzl ±Boomerang;76] ± 2dFDR 6^] data. 

finite horizon at the magnetic field generation time). Based upon this, the fact our results favor a negative spectral 
index near -3, indicate that an inflation mechanism 0-01 niay be the most likely origin for the PMF. 

A PMF affects not only the temperature fluctuations, but also the polarization of the CMB. Although we fit all 
available polarization data, it turns out that the TT and BB modes (where T is the temperature fluctuation and 
B is the curl- like component of polarization) are the most important. Figures [T2] and [T3l show a comparison of the 
computed best-fit total power spectrum with the observed CMB spectrum. Plots show various spectra for the TT and 
BB modes. We plot the best fit and allowed regions both including the SZ effect (scattering from re-ionized electrons) 
at the K(22.8GHz) band (upper curves) and without the SZ effect (lower curves) in FigurdTT] Including the SZ effect 
only slightly diminishes the best fit magnitude of the PMF. Although the CBI point falls about 1 a above the best 
fit, the is dominated by the better ACBAR08 data and this point does not significantly affect the deduced PMF 
parameters. 

DISCUSSION 

Effects of a PMF on the early universe have been studied by a number of researchers including ourselves. A summary 
of some of the relevant topics and their implications for different PMF field strengths is given in Table 2. For example, 
in the first row, if a PMF originated during a post inflation epoch, its strength is probably |Ba| InG, while if 
it was generated during inflation its field strength at the surface of last photon scattering was likely |Ba| ~ InG. 
Furthermore the BB mode of the CMB from the tensor component is strongly affected by the gravitational wave 
background. Therefore, the BB mode also shows characteristics of the gravitational wave background from the PMF 
and an inflation origin. If the PMF has a field strength order of 1 riG, then from figure [5] and [6l we expect that the 
gravity wave background (without the PMF) dominates the BB mode for smaller t < 200, while, the PMF dominates 
this mode for i > 200. Hence, if the gravity wave background without the PMF is generated by inflation, the origin 
of the gravity wave background on large scales will be inflation, while the background for smaller scales is from the 
PMF. 

Based upon the current status of the field, we can summarize the following: 1) A PMF mainly affects the CMB 
temperature fluctuations and polarization anisotropics for small angular scales; 2) Fits to observations of the CMB 
polarization anisotropics for small angular scales are generally better when effects of the PMF are included; and 3) A 
PMF can either increase or decrease the matter density fluctuations on scales less than that of galaxy clusters. 

Comparison with other determinations 

As noted above, in our previous numerical estimate [29^1, the curvature perturbations of the scalar mode were too 
small to stabilize the numerical calculations for large scales and early times. This is caused by an unwanted cancellation 
between contributions from the primordial magnetic field and the primordial radiation fields which seemed to imply 
excess power for the lowest mulitpoles of the CMB. In order to solve this problem, we have adopted the numerical 
approach and initial conditions for the matter contributions from Ref. [79|. This formulation has been developed 
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independently and gives consistent results for large scales and early times. We have now confirmed that our results 
are consistent not only with results of Ref. (79j , but also those of other groups which employ semi- numerical methods 
[e.g. Ref. [3fl]]. 



Influence of the PMF on Large Scale Structure 

In Ref. 0] we deduced the effects of a PMF on the matter energy density fields by considering a stochastic 
PMF that depends upon scale. We then quantitatively discussed the effect of a PMF on the seeds of LSS in the 
early universe. We have also considered more general effects of the PMF than those considered in previous work. 
For example, we included not only the magnetic field tension but also the increases in pressure and energy density 
perturbations from the field. Furthermore, by considering the correlation between the PMF and the matter density 
fluctuations, and utilizing mathematically exact stochastic PMF power spectrum sets, we have obtained reasonable 
and accurate models for the evolution of baryons, CDM, and photons, and therefore, the large scale structure. We 
have shown (ssj that the PMF can play a very different role on the evolution of density perturbations depending upon 
how the PMF and matter fluctuations are correlated. We have also considered the fact that after decoupling, the 
CDM is influenced indirectly by the PMF through gravitational interaction. 

We reported in Ref. [29| that a PMF on small angular scales scales (A «lMpc) and a field strength oi B\ ^ 1 
nG [3, 12, 18, 27, soil provides a new interpretation for the excess of CMB anisotropics for high multipoles. We 
showed that if a PMF with such strength was present, it is very likely that it has affected the formation of large 
scale structure. In Ref. [sij it was suggested that in order to avoid false coupling between the baryons and the CDM 
on small scales, one should use independent transfer functions for the baryons and CDM. A PMF would be another 
source of this difference in the transfer function for baryons and the CDM. Since the density perturbations in the 
early universe have evolved to the present LSS, the evolution of the LSS with a PMF becomes different from that 
without a PMF. We have shown that the baryon and CDM energy density perturbations evolve very differently in 
the presence of a PMF; and with the PMF taken into consideration, the evolution of large scale structure becomes 
more complicated. 



Future Observations 



Magnetic fields affect the generation and evolution of a wide variety of cosmological and stellar objects in the 
universe, e.g. galaxy clusters, galaxies, primordial stars, etc., (cf. Table 2). Therefore, it is important to constrain the 
parameters characterizing the PMF in order to study their effects on physical processes in the universe. In this regard, 
forthcoming missions, e.g. QUIET, PolarBear, and Planck are expected to provide much more precise data on the 
CMB temperature and polarization anisotropics. At the same time data on the formation of large scale structure will 
continue to be measured more precisely both by future space-based observations and from continuing ground-based 
observations such as the SDSS and 2dF projects. 

As summarized in this review, at the present time, we now understand the main effects of a PMF on both the 
CMB and the matter power spectrum. The important task remaining in the future will be to apply our new and 
correct numerical methods for calculating all PMF fluctuations of the scalar, vector and tensor modes in the TT and 
BB power spectra (see figures [TTl [12] and [131) to fit the new observations as they emerge. Only through a critical 



Topics 


|Ba| < InG 


|Ba| ~ InG 


the origin of PMF 


post-inflation 


Inflation 


the gravitational wave 
background 


Inflation 


Large scalcFinflation Small scaleFPMF 


the formation of LSS 


necessity of some generation model for 
the magnetic fleld or vorticity 


necessity for proper initial conditions 
considering the magnetic fleld and 
vorticity 


the magnetic field in 
clusters of galaxies 


necessity of some global magnetic fleld 
generation mode -I- amplification model 


PMF origin 


the re-ionization and 
the first objects 




Formation of the flrst stars and the re- 
ionization are stimulated by a PMF. 



TABLE II. Effects of the PMF on each physical processes 
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comparison with such future data can we hope to ultimately clarify the magnitude and origin of the PMF in the early 
universe. 

Although, some degeneracy between the standard cosmological and PMF parameters is expected in the present 
data, we believe that by simultaneously investigating both the CMB anisotropics and matter power spectra one can 
resolve these degeneracies. The ultimate establishment of the existence and properties of a PMF will permit a better 
understanding of its generation and evolution and also provide new insight into the formation of LSS as well as a 
possible new probe of the physics of the early universe. 

This work has been supported in part by Grants-in-Aid for Scientific Research (20105004, 20244035, 21740177, 
20169444) of the Ministry of Education, Culture, Sports, Science and Technology of Japan. This work is also supported 
by the JSPS Core-to-Core Program, International Research Network for Exotic Femto Systems (EFES). Work at UND 
supported in part by the US Department of Energy under research grant DE-FG02-95-ER40934. 



Initial conditions 



In this Appendix, we introduce the initial conditions of compensated magnetic modes with massless neutrinos in 
the synchronous gauge which were derived by Ref. 56] . Subscripts in this appendix indicate the standard cosmological 
species: 7: photon, v. massless neutrino, b: baryon, and CDM: cold dark matter. We define the density closure 
parameters as the ratios of density parameter as VLx — Px/ Per, where pcr is the critical density of the universe and 
the subscript x denotes the various species listed above. We also define ratios of = fi^/Jlr, R^ = Vt^/Vl^^ and 
i?b = i^b/i^m, where fir = + and i7m = f^cDM + f^b- 

scalar mode 

Potential 

h = -iR^^T + -R^-^t"^ +0{t^) (21) 



"-7 ' n '•7 9 

a 2 
Curvature potential 

1. P 



R-y—T 



2 'a 



R^Ry ,2 3^ P 



6(4i?^ + 15) 4 ''a^ 



21 

+ 0{t^) (22) 



Density perturbations 

S^c^M = -?^7 + - UXr' + 0{r^) (23) 



^ + 2i?/r - ( + 3i?^^ ) t2 + 0{t^) (24) 



-7 - -7 ' --7^' V 6 
a V 6 



= -R^ + 2i?/r + ( - iR,^ ] t2 + 0[t^) (25) 



5':' - (26) 



4 

Velocity perturbations 



i;(s) ^ -R^kT + i(l - 2R^ - SRbRu + Rl)-kT^ + 0[t^) (27) 
'44 a 

1 , 1/5, 



= -R fcr + -i?^-fcT2 + e)(r3) (28) 

4 ' 4 'a 

«f=«f (29) 
Anisotropicstress 

Boltzmann hierarchy of neutrinos 
vector mode 



Vector potential 



56i?^ + 210 7(8i?2 + 90R^ + 225) a ' ' 
Velocity perturbations 

Anisotropicstrcss 

Boltzmann hierarchy of neutrinos 



tensor mode 



tensor potential 
Anisotropicstress 



Boltzmann hierarchy of neutrinos 
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FIG. 1. Comparison of CMB temperature fluctuations from the scalar mode with (red thick curve) and without (thin dashed 
curve) terms of C(t^) from the matter contributions to the scale factor. For this plot we have set B\ — 1.0 nG and ub = —2.9. 
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FIG. 5. CMB temperature fluctuations from the PMF. Plots show various ranges for: (a-1) and (b-1) TT(2 < i < 2750), and 
(a-2) and (b-2) TT(400 <i < 2750) with Bx = 3.0 nG and ns = -2.0 (left panel) or -1.0 (right panel). Green dashed, red 
bold, blue dashed-dotted, and black dotted curves show positive, no and negative correlations, respectively. The black dotted 
curve shows the primary spectrum without a PMF. 
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FIG. 6. CMB polarization anisotropies from the PMF. Blue thin and red bold curves show {Bx,nB) ~ (3.0nG, — 2.0) and 
{Bx, Mb) ~ (3.0nG, —1.0). The black dotted curve shows a primary spectrum with weak lensing effects (without the PMF). 




FIG. 7. Curves of constant values for the as in the parameter plane of PMF amplitude B\ vs. spectral index ub for the case 
of no correlation between the PMF and the matter density fluctuations. Red dashed, blue dotted and green bold curves show 
constant values of as = 1-0, 0.9 and 0.8, respectively. 
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FIG. 8. Curves of constant values for the ug in the parameter plane of PMF amplitude Bx vs. spectral index ub for the 
negative(s[Dp] = —1) and positive(s[op] = 1) correlations between the PMF and the matter density fluctuations. Panels (a-1) 
and (b-1) show ub < —1.5 as the PMF pressure dominated region. On the other hand, panels (a-2) and (b-2) show ub > —1.5 
as the PMF tension dominated region. Red dashed, blue dotted, green bold, azure dotted thin, and black dashed thin curves 
show constant values of as = 1.0, 0.9, 0.8, 0.7 and 0.6, respectively. 
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FIG. 9. Probability distributions and contours of la and 2a confidence limits for the standard cosmological parameters as a 
function of the PMF field strength \Bx\ and power law index ub- Orange contours show 1 a{68%) confidence limits and yellow 
contours show 2 a{95%) confidence limits. Blue curves in the left and bottom of the figure show the probability distributions of 
each parameter. Note the existence of a maximum in the probability distributions for \Bx \ and na- The standard cosmological 
parameters do not have a degeneracy with the PMF parameters because they are mainly constrained by the observed CMB 
data for £ < 1000 (up to the 2nd peak) , while the PMF is mainly influenced by the power on smaller angular scales and higher 
muhipoles, £ > 1000. 
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FIG. 10. Probability contours in the plane of B\ vs. ne. Orange contours show 1 a(68%) confidence limits and yellow contours 
show 2 a(95%) confidence limits. 
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FIG. 11. Probability distributions and contours of la and 2a confidence limits for the parameters as, Ho, ^rcion and Age. 
Orange contours show 1 a{68%) confidence limits and yellow contours show 2 a(95%) confidence hmits. Blue curves show 
probability distributions for each parameter. Note that these are not input priors, but output results. 
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FIG. 12. Comparison of the best-fit computed total power spectrum with the observed CMB spectra. Plots show various ranges 
for: (a) TT(2 < t < 2750) and (b) TT(400 < t < 2750) modes. Orange regions are from the best-fit parameter set and allowed 
regions are from constrained parameter set as Table. 1. Red, black, blue, green, orange, and purple dots with error bars show 
WMAP 5yr, ACBAR 08, Boomerang and CBI data, respectively. Since the SZ effect depends upon frequencies of the CMB 
for the TT mode, we plot the best fit and allowed regions in panels (a) and (b) surrounded by the curves with the SZ effect at 
the K(22.8GHz) band (upper curves) and without the SZ effect(lower curves). 
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FIG. 13. Comparison of the best-fit computed total power spectrum with the observed CMB spectra. The orange region is 
from the best-fit parameter set and the allowed region is from the constrained parameter set as given in Table. 1. Red, black, 
blue, green, orange, and purple dots with error bars show WMAP Syr, ACBAR 08, Boomerang, CBI, CAPMAP, and DASI 
data, respectively. Downward arrows for the error bars on this figure indicate that the data points are upper limits. Since 
the BB mode is not affected by the SZ effect, the upper curves and lower curves in this figure are defined by the constrained 
cosmological parameters of Table. 1. 




